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Abstract 

Stationary solutions of the Chern-Simons effective field theory for the frac- 
tional quantum Hall systems with edges are presented for Hall bar, disk and 
annulus. In the infinitely long Hall bar geometry (non compact case), the 
charge density is shown to be monotonic inside the sample. In sharp contrast, 
spatial oscillatory modes of charge density are found for the two circular ge- 
ometries, which indicate that in systems with compact geometry, charge and 
current exist also far from the edges. 
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I. INTRODUCTION 



The present work focuses on charge and current distributions in clean two dimensional 
electronic systems with edges which are subject to strong perpendicular magnetic fields. 
Investigating the physics of the integer or fractional quantum Hall effect (QHE), and in 
particular, elucidation of the precise charge and current profiles in these systems is a fun- 
damental problem from both theoretical and practical experimental points of view. The 
quantum mechanical dynamics of electrons in two dimensional systems at strong magnetic 
field is characterized by two fundamental concepts. The first one is the formation of Landau 
levels which play an essential role in the studies of integer QHE. The second, and probably 
more profound is the effect of the Coulomb interaction and the emergence of the fractional 
QHE. 

Indeed, a substantial theoretical effort has been devoted to compute the charge and cur- 
rent distributions within a realistic quantum mechanical picture. For the integer QHE in an 
infinitely long Hall bar, a self consistent formalism relating charge density and electrostatic 
potential has been suggested 0. Numerical solutions of the pertinent equations indicate 
that the charge density is monotonic decreasing inside the Hall bar, and has a power law 
singularity (\x — L\~ l l 2 ) at the edges In the particular geometry of a semi-infinite plan, 
an analytic solution using the Wiener- Hopf method is obtained M, and the charge density 
is found to have the same power law singularity near the single edge. 

For two dimensional electronic systems at very strong magnetic fields the electron- 
electron interaction plays a fundamental role, as it leads to the fractional quantum Hall 
effect (and in fact, according to the global picture |§, also to the integer quantum Hall 
effect). Hence, an evaluation of the charge and current profiles in the fractional QHE is 
somewhat more intricate. One of the successful methods for describing the fractional QHE 
is the effective Chern-Simons gauge field theory which includes the Coulomb interaction 
in a non-perturbative manner. Employing this theory for the calculations of charge and 
current distributions is therefore distinct from various other approximations used for treat- 
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ing these quantities in interacting systems. These include Hartree ||, Hartree-Fock and 
others 0]. 

In the present work, calculation of charge and current profiles are based on the Chern- 
Simons field theory for treating the fractional quantum Hall liquids. It turns out however, 
that in order to preserve gauge invariance, inclusion of edges within this formalism must 
be done carefully. One possible approach || is to add a one dimensional edge term to the 
original action. (See also PJiOj and refs. therein.) Once it is assumed that all the important 
physical quantities are concentrated near the edges, this chiral edge action becomes a pow- 
erful tool for studying the various edge effects. Thereby, the theory is formulated through 
the Tomonaga-Luttinger liquid description, whose experimental verification is under intense 
investigation especially in tunneling experiments. Recently, such tunneling effects were mea- 
sured and found to be consistent with theoretical predictions ||1 1|| . 



While studying the physics pertaining to stationary quantum Hall states in which there 
are 'non-zero' currents flowing, one has no apriori knowledge of whether or not charge and 
current are concentrated near the edge. In order to approach this problem, we start our 
investigation from the Chern-Simons effective theory which treats edge and bulk properties 
on an equal footing |T2||. It leads, among other results, to a set of self-consistent equations 



for charge and current distributions within the two-dimensional sample. In the most general 
case, for filling factor v = m/ (mp + 1) with integer m and even p, there are 3m x 3m coupled 
integro-differential equations. Here we restrict ourself to stationary solutions with m = 1 
(non-hierarchial case). For clean systems with sufficiently high symmetry it ends up with 
a single homogeneous eigenvalue problem determined by an integro-differential equation. 
This equation was derived in ref. |12| for the Hall bar geometry, and appears to be markedly 
different from its anlalogous one for the integer QHE [|TJ . 

One of the interesting problems is to study the spectrum of the eigenvalue equation thus 
obtained. For the integer QHE, it is suggested that the spectrum is continuous [jl]. On the 
other hand, nothing is known for the fractional case so far. If we assume that the spectrum 
is continuous also for the fractional case, the integro-differential equations can be solved 



analytically order by order and it leads to a charge distribution which is finite at the edges 
0- 

We present below a general method for solving the pertinent equations both for the inte- 
ger |I] and the fractional |nj QHE systems. Unlike the previous order by order treatment, 
our method is based on the Green function method for solving eigenvalue problems. The 
general method is then applied for some specific systems of physical interest. 

Our results indicate that for the infinitely long Hall bar geometry, only a single eigenvalue 
exists, which is physically acceptable, and the corresponding charge distribution agrees with 
the one obtained earlier using the analytic method ||13|| . For the disk and annulus geometries, 



many physically acceptable eigenvalues are present, and the corresponding solutions behave 
as radial modes p„(r) with n nodes between Ri and i?2- This surprising appearance of 
spatial charge oscillations in fractional quantum Hall systems with compact geometries has 
not been noticed before, and might shed a new light on the physics of edge channels. 

In Section 2, we explain the Green function method which is used in solving the equa- 
tions for the charge density. For convenience, the method is explained within the Hall bar 
geometry. The solutions for the charge and current profiles in a disk and an annulus, which 
comprise the novel part of the present work are presented in Section 3. A few specific topics 
are discussed in the Appendices. In Appendix A, the Green function method for the integer 
QHE system || is discussed, while in Appendix B, the results for the fractional QHE in the 
Hall bar geometry are presented and shown to be consistent with those obtained within an 
analytic treatment fT3| . Finally, in Appendix C, the Green function method is combined 
with a power expansion technique which proves to be useful for the Hall bar geometry. 

II. THE GREEN FUNCTION SOLUTION 

Consider an infinitely long Hall bar stretched along the y axis between the points x = L\ 
and x = L2 subject to a strong magnetic field H in the z direction. In the stationary QHE 
the charge density p depends only on x, and at present we are interested in the density profile 
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p(x) for non-hierarchial filling fraction v = l/(p + 1) with p even. Here p is the difference 
between the charge distribution and its average. The equations derived in Ref. [|T^] couple 
the charge density p(x) and the current density J y (x), but if the later is eliminated, a single 
integro-differential equation for p(x) remains, which reads 

(S 7 ,ucg- 1 ) 2 -^p(x)-p(x) = (4z/ 2 ^" 1 )^ J2 log\x - x'\p{x')dx' . (1) 

In this equation c is the velocity of the Chern-Simons gauge field and g~ l is the coupling 
constant of the Maxwell term in the Chern-Simons Lagrangian. The constant £ > comes 
from the solution of the Poisson equation relating the electrostatic potential to the charge 
density. Eq.(P must also be accompanied by the condition 

p(x)dx = 0. (2) 

It is useful to define q = (Sirucg^ 1 )^ 1 and scale coordinates x — > qx so that the integration 
limits are Xi = gLj. Eq.([l|) then reads, 

d? d 2 

— p{x) - p{x) = p— log\x - x'\p{x')dx', (3) 

where all the quantities appearing in the above equation are dimensionless, including the 
constant p = v^jl'nc. 

Before proceeding with the solution of Eq.(^) three general remarks are useful at this 
point. 1) Eq.(^) (and similar equations which are obtained later for other geometries) is a 
peculiar eigenvalue problem for the charge density p, which can formally be written as 

Ap + pBp = 0, (4) 

where A and B are certain linear operators and p is an eigenvalue. In Eq.(|j), the operator A 
is self adjoint but B is not. However, physics requires the existence of real eigenvalues. We 
show below that Eq.(|3|) can indeed be cast into an equation having the same form as Eq.(|J) 
in which A is symmetric and B is symmetric and positive definite. This problem has then 
only real eigenvalues. 2) The value of p appearing therein cannot be chosen arbitrarily, but 



must be selected from the relevant set of eigenvalues. On the other hand, following Eq.(^), 
the constant p is composed of certain physical quantities such as mass, charge, dielectric 
constant, filling factor etc. If the spectrum is discrete, it might be tempting to formulate 
here a kind of quantization rule. This is of course too ambitious, since the basic theory is 
not an exact one. On the other hand, the sign of p should be consistent with its physical 
content. 3) It is not useful to perform the second derivative on the right hand side of Eq.([3|) 
and obtain an integral equation. Indeed, the resulting equation will have a very singular 
kernel which turn its solution practically impossible. Instead, we use the Green function 
method which is so successful in solving eigenvalue problems of the Sturm-Lioville type. 

Without loss of generality we can set the Hall bar symmetrically between the points 
X\ = —qL and X2 = qL (namely, L2 = —L\ = L). Furthermore, following Refs. P|JT3]| 
we limit our set of solutions to be antisymmetric, p(—x) = —p(x) which automatically 
satisfy the condition (||). The coordinate x is then limited within the interval [0, qL] and 
the equation for the charge density becomes, 



d 2 d 2 
—p(x) - P (x) = p— 



pqL 


x — x' 


/ log 


/o 


x + x' 



p(x')dx', (5) 



together with the condition p(0) = 0. The value of p(qL) is not specified. 
We look for a (symmetric) Green function G(x, x') = G(x',x) which, formally can be con- 
sidered as the inverse of the operator — 1] with the appropriate boundary conditions. 
Hence, it should satisfy 

d 2 1 

-— — 1 G(x, x') = 8{x — x ), (6) 
dx l 

G{0,x') = 0. (7) 

Denoting respectively by x< and x > the smaller and larger values of x and x', it is easily 
verified that 

G(x, x) = — sinh(x<) cosh(x>) + a sinh(x) sinh(x'), (8) 

where the second term on the right hand side reflects the freedom resulting from the absence 
of a second boundary condition at x — qL. 
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We now express p on the left hand side of Eq.([|) in terms of G(x, x'). Formally we apply 
the operator G = [j^ — on both side of @ and use the formal identity 



dx 2 



1 1_ 

dx 2 



d?_ 

dx 2 



(9) 



Strictly speaking, the operator [ J-^ — 1] has a zero eigenvalue corresponding to the function 
sinha; so that its inverse must be defined within the subspace that is orthogonal to this 
function. In so doing we might abandon some of the solutions of the original equation (^|). 
We will see however that the set of solutions is rich enough to capture the main physical 
content. With this point kept in mind the integro-differential equation (|5|) is transformed 
into an integral equation, 



qL 



p(x) = p I dx log 



x + x' 



p(x') + p dx" dx' G(x, x") log 

J JO 



x" + x' 



p(x'). (10) 



Note that the right hand side of Eq. (|10|) vanishes at x = as required. The kernel of the 
integral equation does not appear to be symmetric, so, apriori, there is no guarantee that 
the eigenvalues p are real. In order to proceed and actually show it, it is convenient to write 
Eq.fllOD in its operator form, 



P = p(l + G)Lp, 



(11) 



where the integral operators have their obvious coordinate representations (x\G\x') 
G(x,x'), (x\L\x') = log 



(not to be confused with the system length) and = 
5(x — x'), all of them are symmetric. It is not difficult to show that the operator — L is 
positive definite, and therefore its square root (— L) 1 / 2 exists. Eq . ( |TT1) is then equivalent to 
the following equation for 7] = (— L) 1 / 2 p, 



V = p(-L)^[-{l + G)](-Lf/\ 



(12) 



whose kernel is symmetric. Only solutions corresponding to positive eigenvalues p are phys- 
ically acceptable. As we have commented after Eq.(|9]), every solution of Eq.(|T2"D is also a 
solution of Eq.(|5|), but the converse is not necessarily true. Numerical solutions of Eq. (|T2"D 
will be discussed in appendix B. 



III. OSCILLATING CHARGE AND CURRENT DISTRIBUTION IN DISK AND 

ANNULUS 



We now use the Green function method introduced above to solve the equations of Ref. 
in a circular geometry and demonstrate the occurrence of spatial charge oscillations. The 
algorithm is demonstrated for a disk geometry, and then, later on, some minor modifications 
are introduced in order to study the annulus geometry We believe that for both systems 
our results are experimentally relevant. 

Consider a clean disk of radius R subject to a strong perpendicular magnetic field such that 
the filling factor is v = l/(p + 1) with p an even number. For systems with axial symmetry 
the charge density p, the electrostatic potential V and the (tangential) current density Jg 
depend only on the radial coordinate r. There is of course no radial current. It is useful to 
define the radial differential operator = + - j-. Then, we obtain the following set of 
equations, 

1 4z/ 2 

P(r) ~ -Dlp{r) = —DlV{r) (13) 
9 9 

V(r) = -f [ (f)(r, r')p(r')r'dr' (14) 
Jo 

0(r ' r )= /n / 2. 3 9 / W (15) 

Jo v + t — *rr cos t) 
where q = (Rirvcg^ 1 )^ . Of course, the last two equations are just the solution of the Poisson 
equation expressing the electrostatic potential V(r) in terms of the charge density p(r) for 
systems with axial symmetry, with £ > 0. The charge density and the electrostatic potential 
are expected to be finite in the disk. In particular they should be regular at the origin. The 
weak singularity of 0(r, r') at r = r' is then integrable. 

Eliminating V(r) from Eqs.(^) and transforming to dimensionless coordinates x = qr we 
get a single integro-differential equation for p(x), 

[Dl - l\p(x) = nDl / <f)(x,x')p(x')x'dx', (16) 
Jo 

where p, = v^flixc > 0. 



We look for a (symmetric) Green function G(x,x') = G(x',x) which, formally can be con- 
sidered as the inverse of the operator \D 2 X — 1] with the appropriate boundary conditions 
(regularity at x — 0). Hence, it should satisfy 

[Dl-l]G(x,x>) = 5 -^±, (17) 

x 

G(0,x') : finite. (18) 

Denoting respectively by x< and x> the smaller and larger values of x and x', it is easily 
verified that 

G(x,x') = -/ (x<)K (x > ) + q ; /o(x)/o(xO, (19) 

where Jo and Kq are the modified Bessel functions. The second term on the right hand side 
reflects the freedom resulting from the absence of a second boundary condition at x = qR. 
We can now repeat the same procedure which led from Eq.(||) to Eq. flToD , and later to 
Eq . (|l2l) , with some slight modification due to the presence of the volume element x dx in 
the relevant integrals. Thus we define the operators G and $ such that their configuration 
space representations are {x\G\x') = yfxG(x, x')y/a^ and = y/x<f)(x, x^y/a?. It is not 

difficult to show that the operator $ is positive definite, so its square root is well defined. 
The function r] = $ 1 / 2 x 1 / 2 p then satisfies the integral equation 

r] = /i$ 1/2 (l + G)$ 1/2 r/, (20) 

whose kernel is symmetric. 

For an annulus with radii R\ < R2 the Green function has much more freedom because the 
origin is not reached. Instead of ( |T9"D we may now have 

G(x,x') = -[a/ (^<)i^o(^>) + ^Kq(x < )Iq(x > )} + jI (x)I (x') + SK (x)K (x'), (21) 

with a + P = 1 . 



We have solved Eqs.(^Tj) numerically for the disk and the annulus geometry using iV = 
160 Gaussian integration points for the radial integration. As a credibility check of our 
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numerical procedure we have confirmed that the numerical value of the positive eigenvalues 
is independent of N . 

The operator appearing on the right hand side of equation (p0| ) is real and symmetric. 
One can then look upon the eigenvalue equation as an equation for a rope which is not 
attached at its ends. This analogy is somewhat expected, since equation QT5D is very similar 
to an inhomogeneous wave equation for stationary solutions. The resulting solutions p n {x) 
corresponding to positive eigenvalues X n = l//x n , n = 0, 1, 2, .. with \ n < X n +i appear to have 
n radial nodes. The lowest mode n = cannot satisfy the condition of zero total charge @, 
and hence it is physically unacceptable. All the other modes indeed have zero total charge, 
and can therefore be regarded as representing the pertinent radial charge density. 
For the disk geometry we chose qR = 1 in Eq.([T6|) and a = in the expression ([19]) for 
the Green function. The charge density profile for the lowest four modes n = 1,2,3,4 is 
displayed in Fig. [I]. The basic features of these solutions are as follows: 1) The charge 
density is maximal at the center. 2) The mode n = 1 with one zero reminds us of the charge 
distribution in the Hall bar geometry although it does not have a sharp peak at the edge. 3) 
Most notably, there are numerous oscillatory solutions which, to the best of our knowledge, 
were not predicted so far. 

For the annulus geometry we chose qR\ = 0.5 and gi?2 = 1.0 as integration limits instead of 
and qR in Eq. flllf ). We also take a = l,/3 = 7 = 5 = 0in the expression (|2l|) for the Green 
function. The charge density profile for the lowest four modes n = 1,2,3,4 is displayed in 
Fig. |2|. The basic features of these solutions are different from the ones found for the disk 
in that they do not have a maximum at the inner edge. From that point of view, a disk 
cannot be considered as an annulus with vanishingly small inner radius. The large value 
of the charge at R = is attributed to the weak (but integrable) singularity of K (x) near 
x = 0. However, here again there are numerous oscillatory solutions. We know and 
shall also see in Appendices B and C that in the Hall bar geometry, there is only a single 
solution, which is monotonic. The novel result of spatial oscillatory charge density must 
then be related to the fact that the pertinent physical systems are compact. It then belongs 
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to the realm of mesoscopic and quantum dot physics, which, nowadays, is experimentally 
accessible. 

We now present our results for the current profile in the circular geometry. The derivation 
of the current once the charge density is given has been explained for the Hall bar geometry 
in refs. [0,|T3J. Here we are interested mainly in the oscillatory pattern. It is not difficult 



to show that the tangential current Jg(x = qr) is given by, 

J$(x) = a— I x p(x')dx + b, (22) 

X J qRi 

where a and b are constants and R\ — for the disk and finite for the annulus. In the 
following we display the integral appearing on the RHS of ( ^2]) since it gives the main 
characteristics of the actual current. In figure ^ the current profile for the disk is shown, for 
the same configuration as in figure |l|. The current of the lowest mode is similar in shape 
to the current evaluated in the Hall bar geometry [|1 ,13|| . The current of other modes is 



oscillating, which shows clearly that for these modes, current is not concentrated solely near 
the edges. Similar result occurs also for the annulus geometry, for which the current of the 
four modes pertinent to figure |2] is displayed in figure |j. 

One of us (Y.A) would like to thank the Japanese Society for the Promotion of Science for 
supporting his visit to the Institute for Solid State Physics in which this work has been 
carried out. 



APPENDIX A 

Solution of the equation for the integer QUE. 

In this appendix we present a Green function solution of the equations derived in Ref. ]l| for 
the integer quantum Hall effect in a Hall bar. If the Hall bar stretches along the y direction 
between points x = L\ and x = L 2 , the charge density and the electrostatic potential depend 
only on x. They are determined by a set of equations, 
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d 2 fi L * 

p{x) = / log 
dx z Jo 



p(x')dx', (25) 



V{x) = -2e 2 f L2 log \x - x'\p{x')dx' (23) 
p(x) = j^-V"(x), (24) 

where n is the Landau level number and uo c is the cyclotron frequency. Transforming to 
dimensionless coordinates x — > qx this set of equations is equivalent to a single equation for 
p(x) which, for antisymmetric solutions reads 

d 2 f qLx 

Ioet 

x + x' 

where p = (4e 2 nq)/(huj c ) > and L x = L2 — —L\ is half the width of the Hall bar. Here q 
is just a parameter with dimensions of inverse length. In order to apply our Green function 
method for solving Eq.(^) we add p" on both sides and use algebraic manipulations as in 
Section 2. The result, using the notation of section 2, is 

-Gp = p(l + G)Lp. (26) 

If we apply the operator (1 + G)^ 1 on both sides, we get 

[l-(l + G )- 1 ] P = p(-L)p. (27) 



Equation (p7|) is of the form (P with symmetric A and symmetric positive definite B and 
therefore it has real eigenvalues p. It would have been more useful to solve an equation for 
the potential V(x) rather than the density p(x), since the later one is known to have a power 
singularity near the edges [p],Pi|13l. Unfortunately we were unable to obtain an equation of 
the form @ with symmetric A and symmetric positive definite B. 

The charge density in a Hall bar for the integer QHE is displayed in Fig. |5|. For the actual 
numerical solution we take the Hall bar to be located between x = — 2 and x = +2 in order 
to compare the charge density of the non-interacting system with that of the interacting 
system reported below in Appendix B. As can convincingly be inferred from Fig. [| the 
charge density in the integer quantum Hall regime is singular near the edge. This is of 
course consistent with earlier results |Tl,p|JT3[| . 
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APPENDIX B 



Charge density profile in the Hall bar geometry for the fractional QHE. 
As was stated already, the charge distribution profile in the Hall bar geometry for the 
fractional quantum Hall system was calculated by an analytic iterative method ||13|| . In 
this appendix we present the exact solution obtained within the Green function technique 
introduced in Section 2. We have solved Eq.(|T2"D numerically using mesh of N = 400 
integration points between and 2, replacing the operators by N x N matrices. The weak 
logarithmic singularity at x — x' is of course integrable. Out of the N eigenvalues only one 
is found to be positive. This feature as well as the numerical value of the positive eigenvalue 
persist independent on the number N (provided it is sufficiently large). The parameters qL 
(upper limit of integration in Eq.(|T0"D and a (the free parameter in the Green function @ 
should be chosen in such a way that Eq.(|I~2"D has at least one positive eigenvalue. In all 
our attempts, we could not obtain more than a single positive eigenvalue, a result which 
is markedly distinct from that obtained in the circular geometry. For the special choice 
qL = 2 and a = the charge density in a Hall bar is displayed in Fig. || It is reassuring to 



note that this solution is extremely close to the one obtained by iteration methods [|13|] or 
the one obtained by solving Eq. (|i~2"D using a power basis technique as explained in the next 
Appendix. 



APPENDIX C 

Solution of the equations in the Hall bar using a power series expansion. 

In this appendix we present an alternative formulation to study the integro-differential Eq. (^) 

in the form 

u(x) — au"(x) = —/i J dx'ln \x — x'\u"(x'), (28) 

with the antisymmetry condition u(—x) = —u(x). Here, we have used the notation p(x) = 
u"{x). We can then choose {x 2k+1 \k = 0, 1, • • •} as an expansion basis. The integration 
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operator on the RHS of (pq) is represented as 

L ■ x 2l+1 = - ^ dx' In \x - x'\x' 2l+1 

OO c\ 

= y 

^ Q (2k + l){2k-2l-l) 
whose matrix representation reads 

/ n 111 \ 



x 



2fc+l 



(29) 
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2 
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15 
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2 


15 


5 


5 


15 


2 


2 


2 


2 


35 


21 


7 


7 



V 

and the double derivative is 

D 2 x 2l+l - 



(30) 



d 2 
dx 2 



^—x 2l+1 = {2l + l){2l)x 2l -\ 



Thus we have 

(1 - aD 2 )u = fiLD 2 u. 

Since the double derivative D 2 has rapidly increasing next-diagonal elements (21 + 1)(2Z), it 
seems impossible to simulate the eigenvalue problem (|28| ) numerically with any finite size 
truncation of D 2 . A remedy for this is to integrate twice the original equation (|28|) . Let us 
introduce the double integration operator I 2 by 

1 



2^2i+l 



Vx 



CiX + 



-x 



21+3 



(2/ + 3)(2/ + 2)' 

where q is an integration constant. Then all the solutions of equation 

I 2 ^ 2 ^- 1 ^ 2 -a)u = nu, 

also satisfy equation (|28|). After finding some suitable set of constants q, we obtain a positive 
eigenvalue [i whose eigenvector is consistent with the one obtained by the Green function 



method and by an analytic argument |13 
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FIGURES 

FIG. 1. Charge distribution in a disk. The parameters are qR = 1 and a = (see equation 
fll9|)). The four modes correspond to the four positive eigenvalues A„, = l/fi n of equation (|20|). 
with n = 1, 2, 3, 4 which also counts the number of radial nodes. 

FIG. 2. Charge distribution in an annulus. The parameters are qR\ = 0.5, qR2 = 1 and 
a = l,/3 = 7 = <5 = (see equation (|2l|) ). The four modes are as in figure |l[ 

FIG. 3. Current distribution in a disk. The parameters are qR = 1 and a = (see equation 
(|l9|)). The four modes correspond to the four positive eigenvalues \ n = l/fi n of equation (j20|) , 
with n = 1, 2, 3, 4 which also counts the number of radial nodes. 

FIG. 4. Current distribution in an annulus. The parameters are qR\ = 0.5, qR2 = 1 and 
a = l,/3 = 7 = <5 = (see equation (|2l|)). The four modes are as in figure |3[ 

FIG. 5. Charge distribution in a Hall bar for the integer quantum Hall system obtained 
by solving equation (|27|). The Hall bar is located between x = —2 and x = 2 and the charge 
distribution is antisymmetric. The parameter a in the Green function (Q) is set equal to zero. 

FIG. 6. Charge distribution in a Hall bar for the fractional quantum Hall system obtained 
by solving equation (|i~2|). The Hall bar is located between x = —2 and x = 2 and the charge 
distribution is antisymmetric. The parameter a in the Green function (Q) is set equal to zero. 
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